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Abstract 

We propose a technique that allows to simultaneously perform universal control of the evolution 
operator and compensate for the first order contribution of an arbitrary Hermitian constant noise. 
We show that, at least, a three- valued Hamiltonian is needed in order to protect the system against 
any such noise. This technique is illystrated by an explicit algorithm for a control sequence that is 
applied to numerically design a safe two-qubit gate. 
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I. INTRODUCTION 



Within the last decades, quantum control has emerged as one the most fruitful fields in 
both theoretical and experimental physics . In particular, it is of crucial importance 
in quantum computation Q|. To process the information stored in the computer state, one 
must indeed be able to generate any prescribed unitary evolution operator on the computer 
— or, at least a universal set of such operators, unaffected by quantum errors arising from 
the interaction with the environment. 

Different strategies have been developed to control the evolution of closed quantum 
systems, including optimal control approaches and algebraic methods 
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3. 



To deal with open quantum systems, schemes have been designed to cor- 



rect/avoid the undesired effects due to the environment. Quantum error-correcting codes 



(QEC) 
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15! . Q and approaches based on the quantum Zeno effect 
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redundancy of encoding as a way to recover information after the errors occur. Topological 
protection jl9] takes advantage of the symmetries of the system to safely store informa- 
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22|. An alternative to QEC that is 
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25|]. In DD one 



tion in so-called Decoherence Free Subspaces 120, 
substantially less resource-intensive is dynamical decoupling (DD) 
applies a succession of short and strong pulses to the system, designed to stroboscopically 
decouple it from the environment. Similar in spirit to DD, but more general, is the method 
we term here "dynamical control by modulation" (DCM), wherein one may apply to the 
system a sequence of arbitrarily-shaped pulses whose duration may vary anywhere from the 
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In 



stroboscopic limit to that of continuous dynamical modulation 
the DCM approach, the decoherence rate is governed by a universal expression, in the form 
of an overlap between the bath-response and modulation spectra. In such methods, it is 
however not clear a priori whether one can, at the same time, protect the system from noise 
and perform universal control of its evolution operator. 

In this paper we investigate this question by asking whether it is possible to perform 
any arbitrarily chosen evolution of a quantum system while compensating for all Hermitian 
static noises. To be more specific, our goal here is to show how to design a time-dependent 
control Hamiltonian which can, at the same time, impose a chosen evolution to the system 
and eliminate the first order action of any Hermitian static noise. We show that, contrary 
to strict evolution control problems, this objective cannot be achieved with only two- valued 
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Hamiltonians but requires the use of at least three control operators. We go on to show 



le system to evolve under the action of noise only, 



that even a null third operator, causing t 
is enough. Inspired by a previous result [3J|, we propose a new algorithm able to compute 
the appropriate control sequence for any given desired evolution. As an application, this 
algorithm was run to design a safe CNOT gate. 

The paper is structured as follows. We first set the problem to solve and give the ex- 
plicit conditions the evolution matrix must fulfill. Then we show that these conditions are 
impossible to meet by a two-valued control Hamiltonian. Slightly modifying the two-stage 
procedure by allowing for extra steps during which the control Hamiltonian is set to zero, 
we describe an algorithm able to compute an appropriate protected control sequence. An 
application to a two-qubit gate is finally proposed. 

II. CONDITIONS FOR FIRST-ORDER NOISE ELIMINATION 

Let us consider an iV-level quantum system whose Hamiltonian consists of a control- 
lable part, denoted by H c (t), and an unknown but static noise contribution, which can 
be written as a linear combination Af = he^Gi of the Hermitian traceless generators 
{Gi, i = 1, . . . , N 2 — 1} of su(N). The evolution matrix satisfies the dynamical equation 

•ft^p = (H c (t)+M)U(t) (1) 
U(0) = I. (2) 

Upon transforming to the interaction picture relative to H c , one isolates the evolution U (t) 
due to the noise only: defining 

-t 



U e (t)=Texp^J H c (s)dsj 



(3) 



the evolution induced by the control Hamiltonian alone, where T denotes the chronological 
product, one sets U (t) = U\ (t) U (£), which satisfies 

ih— = [Ul(t)NU c (t)]U(t) (4) 



J2^Ul(t)G t U c (t) 



U (t) (5) 



17(0) = I. (6) 
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The first order contribution of the noise to the evolution is thus given by the second term 
in the Dyson expansion of U (t) for the accumulated action, that is 

(t) = V £l f U\ (s) dU c (s) ds. (7) 
i J o 

Our goal is to design a control Hamiltonian H c (t), such that, at the end of the control 
sequence, say at time T c , the evolution operator takes an arbitrarily prescribed value Ud G 
SU (N) while the first order contribution of any constant noise vanishes. We thus require 

U C (T C ) = U d (8) 

and, for any set of constants {£{}, 

(T c ) = V £i f C XJ\ (s) G t U c (s) ds = 0, (9) 
Jo 



that is 

V«, / Ul{s)GiU c {s)ds = 0. (10) 
Jo 

III. ALTERNATING TWO- VALUED OPERATOR SEQUENCE 

Let us now focus on the two-valued alternating perturbation approach. Namely, the 
control Hamiltonian H c (t) alternates between two values A and B for adjustable timings 
which play the role of control parameters. Formally, H c (t) then assumes the bilinear form 
H c (t) = ha (t) A+h/3 (t) B, where a (£) and (3 (t) are two piecewise constant functions taking 
the values 0, 1 and adding up to 1. In the sudden approximation, the overall evolution 
operator induced by such a K-step control sequence has the pulsed form 

U c (T c ) = e ~ lTKHK x ... x e~ lT2H2 x e~ lTlH \ (11) 

where Hk = A when k is even, B when k is odd, and J2k=i = T c . 

Provided that A, B together with their all-order commutators span su(N) (the bracket 
generation condition), it can be shown [3l| that it is possible to design K oc iV 2 such control 
timings {Ti, . . . , T K } for which eq.((8j) is met, i.e. such that 

e- lT >< H « x ... x e~ lT2A x e~^ B = U d . (12) 
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It turns out, however, that, except for the particular case Ud = I, eq. fllOl) can- 
not be satisfied for all noises. As we shall now show, there indeed always exists a 
D-dimensional uncorrectable subspace (1 < D < N — 1) spanned by the noise operators 
C ne /y = \ [A + B, (B — A) n ]. To prove this, let us first introduce the operators 

Uq = I, XJ K >i>i = e~ lT ^ x ... x e~ lT2A x e~^ B (13) 

which allow us to write 

U c (t) = e-^-^^^+'Ui (14) 

(forEi^<t<Ei +1 ^)- (15) 

Let us further define the superoperators 7ij, whose action is given, on any operator X, by 
HiX = [Hi,X\. We can thus write e +eH 'Xe' £H ' = J2T=oi^i X - Since [A,(B-A) n ] = 
[B, (B - A) n ] = [H h (B - A) n ] for every i, we can set C n = H^B - A) n . 
Let us now evaluate the first order effect of C n according to eq.(|T0l): 

/ ' Ul{s)C n U c {s)ds (16) 
Jo 

V / dsUje lsH ^Hi(B - A) n e- lsH *Ui (17) 
i J ° 

ds^-^W-^iB - AYUi (18) 

■o k=0 

iJ2 U i [e lTlHl (B - A) n e~ %TiHi - (B - A) n ] U { (19) 

i 

i^uUB-Aru^-iJ^uHB-ArU (20) 

i i 

iU\T c )(B-A) n U(T c )-i(B-A) n . (21) 

The first order contribution we have just obtained does not depend on the specific timing 
parameters Tj's but only on the final operation U(T C ); it is moreover nonzero for any Ud ^ 
I. Finally, according to the Cayley-Hamilton theorem, the matrix (B — A) cancels its 
characteristic polynomial, which implies that 1 < dim [Span ({(B — A) n , n G N})] < N — 1. 
As a consequence the dimension D = dim [Span ({C n , n e N})} of the uncorrectable subspace 
satisfies 1 < D < N- 1. 

Let us now examine how the two-operator Hamiltonian scheme can be modified so that 
both conditions eq. (l8lfT0l) are met. Suppose we have a sequence of timings {£]_,...,£#}. Now 



Uq and U n >\ are the evolution operator at the end of the n th step of the control sequence 
as defined in eg. ([TBI , while U c (t) is the evolution operator at time t as defined in eg. ffT5t) . 



Following the method described in |3jJ, the last N 2 — 1 timings {tK-N 2 +2, ■ ■ ■ ,Tk} can be 
chosen such that they achieve the evolution Ud ■ U' K _ N2+1 , thus satisfying the condition 
of eq.(E|). The first order contribution of any noise Gi can be explicitly expressed by the 
accumulated action 

K 

Gi^Y^Ul^nUn, (22) 



n=l 

where g^ n = J^ n e-™ Hn Gie +lsHn ds. 

Let us, at each commutation between A and B, allow for the waiting time Ti=i t ... t K during 
which no perturbation is applied. This amounts to adding a third value C = to the control 
Hamiltonian H c {t). We see that the overall evolution operator remains unchanged, while 
the first order contribution of any noise Gi is added the term Yl n =i T nUlGiU n , which is a 
linear function of the waiting times. 

Let us choose the timings {t\, . . . , tx-N 2 +i} such that the coefficients of the waiting times 
span the entire iV 2 • (iV 2 — 1) space. Generally this can be achieved by randomly choosing 
the timings. Thus, by solving a simple set of N 2 • (N 2 — 1) linear equations of the form 
^i~^~J2n=i T nFi,n = {Ti <n = U^GiUn) to find the waiting times t±, . . . , tk, one can eliminate 
the first order contribution of all the noises Gi added during the A, B control sequences. In 
practice, since the waiting times T{ must be positive, this requires a slightly larger K and 
the use of linear programming methods or other minimization techniques. 



IV. ALGORITHM FOR A FOUR-STATE SYSTEM 

This algorithm was applied to a model four-state system, represented in Fig. Q], which 
can be used to store and process two qubits of information. The four states correspond 
to two different angular momenta I = 0,1: |0) = \l — 0,mi — 0), |1) = \l — l,mi — — 1), 
|2) = \l = l,mi = 0) and |3) = \l=l,mi = 1). 

This system is subject to a resonant electric and a static magnetic fields: the 7r-component 
of the electric field couples |0) to |2) while the cr + --components couple |0) to |1), and |0) 
to |3), respectively; the x, ^-component of the magnetic field couples |1) to |2) and |2) to 
|3) while its z-component shifts |1) and |3) out of resonance. Finally, in the rotating wave 



6 



approximation (RWA) , the total Hamiltonian of the system assumes the form 



H c = 


e_cr_ + e a + e_, 


.0+ + b ± A ± + £> 2 A Z 


(23) 


o — — 


11^ Ml 
I 1 / \ u l 


4- \0) (1\ 


(24) 


cr = 


|2> (0| 


+ |0> (2| 


(25) 


a + = 


13) (0| 


+ |0> (3| 


(26) 


A± = 


|1) (2| + |2) (3| 


+ |2) (1| + |3) (2| 


(27) 


A z = 


|3> (3| 




(28) 



where e-^.f and b± :Z are five independent parameters, proportional to the electric and 
magnetic field amplitudes, respectively. By choosing two different sets {e^ 0+ ,bj_ 2 } and 
{ e - o z ] of such parameters, we can define two values for the Hamiltonian H c = A, B 
which may be used as our alternating perturbations. 

After verifying that A and B satisfy the bracket generation condition, we first calculate the 
K = N 2 = 16 timings jtj^ ^ j which realize the identity matrix, by the method described 
in 31|, based on statistical properties of the roots of the identity. By the Newton iterative 

method, we then compute the K = 16 timings {tk=x,...,K} which perform the transformation 

( I o\ 
10 
1 
\0 1 0/ 

Finally, we look for the K (N 2 — 1) = 240 waiting times Ti=\ ,a-(jv 2 -i) which mini- 
mize all 15 possible noise operators to zero. By applying the complete 480-step control 
sequence, consisting of the (N 2 — 1) = 15 repetitions of the sequence {tk=i,...,K=w} plus the 
K (N 2 — 1) = 240 waiting periods {77} with zero control after each step, we can thus impose 
a safe CNOT gate on the system in the presence of any noise that varies slower than the 
control sequence. 



CNOT^-i = C NOTte, where CNOT 



V. DISCUSSION 

It is plausible that the conditions eq. (l8fT0l can be satisfied by an arbitrary three- valued 
Hamiltonian, with C ^ 0. We have looked for an appropriate sequence of timings with such 
a Hamiltonian by direct optimization of a cycle of KN 2 operations and found satisfactory 
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l ± |l> = |/= 1, m,=-l)/ ^ S |2> = |/=1, m, = 0)/ ^ s \3) = \l=l, m,= l) l 



' A, 



A, 1 



/ 

/ 



1 = 



|0> = |/= 0, m, = 0) 



Figure 1: The two-level atomic system (I = 0,1). The arrows show the different couplings due to 
the resonant electric field (operators <t_ j o,+) and the static magnetic field (operators Aj_ iZ ), 

numerical solutions up to = 16. This approach is however slower than the scheme 
described above. 

There are several open issues regarding the universal control method proposed here. 
First, it is important to treat higher orders in the perturbation expansion of the noise and, 
in particular, explore the feasibility of a control Hamiltonian canceling the second order 
contribution of the noise, i.e. the integrals 

( C dt f dsUt{t)G m U c (t)Ul(s)G n U c {s) (29) 
Jo Jo 

for any (m, n). Second, the issue of time-dependent noises is important. The present method 
holds for noise that slowly varies with time, compared to the control sequence. The approach 
must change altogether if fast noises affect the system. Finally, it is imperative to establish 
whether redundancy can be combined with dynamic control techniques to safely process the 
information, as, for example, in an ensemble of identical systems. 
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To conclude, we have raised the question whether universal control of evolution may be 
performed while compensating for the first order contribution of arbitrary constant Hermi- 
tian noise, by means of an alternating perturbation procedure. We have demonstrated that 
this can not be achieved by a two-valued control Hamiltonian: in that case there always 
exists a subspace of uncorrectable noises. If, however, we allow for waiting times, during 
which the system is only subject to noise, our objective becomes feasible. This has been 
demonstrated by an explicit algorithm that yields the appropriate control sequence, as tested 
on the case N — 4. The total waiting time is of the order of half the control period. More 
generally, we numerically checked that a three- valued Hamiltonian can control the evolution 
operator and protect it against any static Hermitian noise. 
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